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The BN/KL region in the Orion molecular cloudB is an archetype in the 
study of the formation of stars much more massive than the Sun!. This region 
contains luminous young stars and protostars, but it is difficult to study because 
of overlying dust and gas. Our basic expectations are shaped to some extent 
^ ■ by the present theoretical picture of star formation, the cornerstone of which 
Q\ ' is that protostars acrete gas from rotating equatorial disks, and shed angular 
momentum by ejecting gas in bipolar outflows. The main source of the outflow 
in the BN/KL region BoB may be an object known as radio source li, which is 
commonly believed to be surrounded by a rotating disk of molecular material 
§'!!. Here we report high-resolution observations of silicon monoxide (SiO) and 
water maser emission from the gas surrounding source I; we show that within 60 
AU (about the size of the Solar System), the region is dominated by a conical 
bipolar outflow, rather than the expected disk. A slower outflow, close to the 
equatorial plane of the protostellar system, extends to radii of 1,000 AU. 

Q-f 

The BN/KL region radiates an estimated 5-8 xlO 4 ergs -1 @0 and is a rich source 
of molecular line emission!0 embedded within the Orion Molecular Cloud (OMC-1). 
Two molecular outflows originate in BN/KL, relatively close to source I, a low- velocity 
(18 kms -1 ) outflow to the northeast-southwest0'@, and an orthogonal high- velocity 
>< '• (30-100 kms" 1 ) outflow0ffi0, both of which extend over about 1' (0.14 pc). Source I 

has a broad-band radio spectrum consistent with an H II region about 60 AU in diameter^!, 
though it has not yet been specifically detected at infrared wavelengths. Masers are often 
associated with massive young stars, and many are broadly distributed across BN/KL 
(masers are the microwave equivalent of lasers). However, a distinct population of SiO and 
H 2 masers is clustered around source To understand better the physical nature of 
source I and its relationship to the outflows in BN/KL, we have mapped these masers with 
high spectral and angular resolution. We resolve the angular structure of the SiO emitting 
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region in detail at each observed line-of-sight velocity, in contrast to previous observations 
with connected-element radio interferometersl'IHIll. 

Maser sources are valuable tracers in molecule-rich regions because they comprise many 
compact, high surface-brightness, spatially distinct points of radiation (spots), with various 
line-of-sight velocities. The spots are test "particles" tracing the three-dimensional velocity 
field of the material in which they are embedded. Maser amplification requires relatively 
small gradients in line-of-sight velocity to achieve long gain-path lengths. To pump the 
maser levels by collisions, the gas temperature must also be at least a few hundred degrees, 
and the density orders of magnitude greater than that sampled by thermal molecular 
emission, i.e., n# 2 ^ 10 8 cm~ 3 . Maser emission is a selective tracer but in star forming 
regions, it can in principle occur in either outflows or disks0. 

We observed emission from the v — 1, J =1-0 transition of SiO (43.12208 GHz rest 
frequency) toward the BN/KL region on 9 July 1995 with the Very Long Baseline Array 
(VLBA) of the NRAO. The data were processed with the VLB A correlator, which provided 
0.22 kms" 1 wide spectral channels, and reduced with standard techniques!!!. Within 
30 kms -1 of the 5 kms -1 systemic velocity^ of the region near source I, we detected 1,212 
maser spots with flux densities between 0.09 and 300 Jy. 

Spectroscopic images show that the SiO maser emission is confined to four regions that 
can be envisioned as marking the arms of an X, subtending about 0''25 (Figure la). We 
propose that the X marks the intersection of two outflow cones with the plane of the sky and 
separates the low-velocity (northeast-southwest) and high-velocity (northwest-southeast) 
outflows in BN/KL. Source I was not detected by us because its surface brightness is 
too low. However, to estimate the position of source I on our maps, we compared the 
flux-weighted emission centroid for each VLBA spectral channel to the earlier, lower 
resolution (0"25) map of Menten & Reidi who measured the position of source I relative to 
the maser emission centroids with 0''01 uncertainty, using the Very Large Array (VLA). The 
relative positions of the dominant clumps of masers in the two maps agree to 0''01. Source 
I lies at the center of the X, within the measurement errors (Figure la). The innermost SiO 
maser features are separated by about the 60 AU diameter of the source I H II region. 

We have added to our study archived H 2 maser data collected on 25 August 1983 
with the VLA (Figure lb), as well as published maps0 of the H 2 maser distribution on 
5 August 1991. The distribution of the H 2 masers is elongated in the direction of the 
low- velocity outflow in BN/KL. Over eight years, the length of the distribution grew by 
about 0''2. We measured proper motions for 19 individual maser features that persisted 
over the eight year period (Figure lb). The motions are primarily parallel to the axis of the 
low- velocity outflow, and the average transverse motion corresponds to about 20 kms -1 , 
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approximately the speed of the large scale low-velocity outflow in BN/KIH 

The previous widely accepted model, in which the SiO masers lie in a rotating, 
expanding, and inclined disk@oi, does not fit our high-resolution VLBA maps. One 
attraction of this model was that the disk could feasibly extend to larger radii, explaining 
the low- velocity BN/KL outflow. In such a disk, the SiO emission is expected to trace 
dynamic surfaces where the derivative of the line-of-sight velocity is small, so that the 
maser gain is high. For most rotation laws, each arm of the X would be expected to show a 
velocity gradient along its length; but three of four do not. Moreover, the distribution of 
SiO masers that we observe shows approximate mirror symmetry in position and velocity 
about a northwest-southeast axis, which constrains the magnitude of any rotation to be less 
than a few kms -1 . If the disk is gravitationally bound to the protostar, then this rotation 
speed, at 25-60 AU radius, constrains the protostellar mass to be less than 1 M . If radial 
forcing enlarges the disk (e.g., a wind) from an initially bound state (with about 1 AU 
radius), then conservation of angular momentum requres a stellar mass less than a few M Q . 
These relatively small stellar masses are inconsistent with the presence of an H II region. 

We propose that the SiO maser emission outlines the limbs of a bipolar, conical, 
high- velocity flow 25-60 AU from the protostar (i.e., the X), and that the H2O masers lie 
in a low- velocity flow in an equatorial band bisecting the high- velocity flow (Figure lc). 
The position angles of the cones are —40° and 150°, with opening angles of 60° and 80°, 
respectively. For an outflow axis close to the plane of the sky and a finite cone wall 
thickness, the longest maser gain paths naturally lie along the limbs of the cones. The 
10-15 kms -1 difference in mean Doppler velocity between the red- and blue-shifted SiO 
masers may be explained if the cones are tipped out of the plane of the sky. The known 
distribution of H 2 emission in BN/KL§ shows two lobes with line-of-sight velocities Doppler 
shifted by more than 30 kms -1 that subtend the openings of the cones. The northwest lobe 
is equally prominent in red- and blue-shifted emission, which suggests that the inclination 
is small. (Blue-shifted emission from the southeast lobe appears to be blocked by a dense 
foreground cloud§§.) 

The dense gas in which the SiO masers lie is probably not gravitationally bound to 
the protostar because there is no detectable rotation. The SiO maser emission probably 
arises at a steady-state shock interface, perhaps between a faster and a slower component of 
the stellar outflow. Such shocks serve as efficient mechanisms that collisionally pump SiO 
masers and reduce depletion of gas-phase SiO onto dust grains. To quantify the lifetime of 
the visible structure, we compared the present SiO maser maps with a partially complete 
VLBI map obtained in 1985c3. The loci of red-shifted masers have remained fixed with 
respect to source I to about 0"01 over 10.3 years, which corresponds to a proper motion of 
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less than 4.5 AU or 2 kms" 1 . We speculate that collimation of the bipolar outflow begins 
at radii less than 25 AU, perhaps inside the HII region, because source I lies close to the 
apices of the outflow cones. 

We fit the H2O maser positions and Doppler velocities to a constant-velocity outflow 
centered on source I, where the H 2 masers are constrained to lie on a toroidal surface, 
as suggested qualitatively by the curvature traced by several H 2 features (Figure lc). 
This surface may be where large-diameter dust grains begin to form and to be accelerated 
by radiation pressure. (Dust condensation and H2O maser excitation together require 
temperatures of 400-1,000 K. If this material is in thermal equilibrium, then the luminosity 
of source I would be 3> 10 4 L .) The torus has the following parameters: mean outer 
radius 360 AU, thickness 270-320 AU, outflow velocity 13 kms -1 , and rotation less than 
a few kms -1 . The torus is tipped down about 10° with respect to the line of sight, as 
evidenced by the southeasterly position bias of the distribution of blue-shifted H2O masers. 
(The masers lie close to the equator.) We assume that the axes of the high- velocity outflow 
and the torus are aligned (Figure lc), from which we infer that the high- velocity outflow is 
at least 50 kms -1 . This limit is consistent with the observed three-dimensional motions 
of Herbig-Haro objects in and the Doppler velocities of thermal CO and SiO 

emissionEJ. We note that the SiO maser emission from the limbs of the cones suggests 
material flowing away from the observer toward the northwest, opposite that expected given 
inclination of the torus. However, we consider that the masers trace gas selectively. For 
small inclinations, optimal amplification can occur with the opposite sense of Doppler shift 
(i.e., red- vs blue-shift) from gas along the axes of the cones, depending on the radial and 
angular dependence of the outflow velocity. 

We conclude that the equatorial region is dominated by outflowing material and not 
by (rotating) natal material. We may extend this outflow to radii of up to 1,000 AU by 
considering the previously mapped distribution of v=0 SiO emissioniia, to which three 
components contribute: high- velocity thermal, low-velocity thermal, and low-velocity maser 
emission. Both low-velocity components overlie the H 2 maser distribution in velocity and 
sky position within 350 AU of source I. The brightness distribution is asymmetric with 
respect to source I, as is the H2O maser torus. On the smaller side of the torus, stronger 
SiO emission coincides with a richer distribution of H2O masers, which is consistent with 
a density gradient in the surrounding medium, as seen at 8.7 and 12.4 /imi'l. Although 
the apparent mixing of H 2 maser and SiO emission is difficult to understand (since 
silicate-grain formation depletes gas-phase SiO), we note that similar mixing is apparent in 
the envelopes of late-type stars@@. 

We note three propitious alignments of structure that strengthen the case for the 
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proposed model. First, no H 2 masers lie within the high-velocity flow demarcated by 
the SiO masers (Figure lc). This flow is probably too hot or turbulent to support maser 
action. Second, the four components of IRc20 shown in Figure 2 subtend the opening 
angle of the northern outflow cone. They may have a direct line of sight to the protostar 
and may be heated by it or by interaction with the high-velocity flow. Infrared light from 
the southeastern outflow is blocked by extinction in a foreground cloud core^S. Third, 
the magnetic field direction inferred from dust emission at 100 /im and 3 mm wavelength 
is parallel to the high-velocity outflow axis@S This suggests that the formation of the 
protostar and the alignment of the outflow axis may have been affected by the broad field 
distribution in OMC-1. 

We argue that the masers clustered around source I provide evidence for a two- 
component outflow from a young star. In the context of star formation theory for high-mass 
stars, no current model readily explains wide-angle collimation of a high-velocity bipolar 
outflow within tens of AU of a protostar and a simultaneous low-velocity equatorial outflow. 
Magnetohydrodynamic (MHD) mechanisms such as X-winds0 and disk winds0 may be 
reasonable candidates. Alternatively, a spherical stellar wind that mixes with and ablates 
an (unseen) accretion disk up to about 350 AU in radius might generate a slow, dense 
equatorial wind, but this would not explain the observed conical structures. In addition, 
the presence of an H II region may preclude a disk. If accretion and outflow processes are 
linked, then the birth of an HII region could dissipate the observed equatorial outflow in 
about 100-200 years. Based on the alignment of the outflows from source I with other 
structures in BN/KL, we suggest that source I drives a substantial portion of the bulk 
motion in BN/KL, though other sources, such as n (Figure 2), rightly draw attention!'!. 
We anticipate that VLBA maps of SiO and H 2 maser proper motions and SiO maser 
polarizations will clarify the structure, dynamics, and magnetic field geometry in this 
region, which will help resolve the debate about sources of energy and outflow in BN/KL. 
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CAPTIONS 

Figure la The v = 1, J = 1-0 SiO maser distribution in the BN/KL region, obtained 
with the VLBA. Unlike previous maps made with connected-element radio interferometers, 
this one resolves the angular structure of the SiO emitting region at each observed 
line-of-sight velocity. The cross represents the centroid and position uncertainty for 
continuum source I with respect to the masers. The double arrows indicate the direction 
of the proposed high-velocity outflow from source I, and the approximate direction of the 
high-velocity outflow in BN/KL. Symbol color indicates Doppler velocity in kms -1 with 
respect to the Local Standard of Rest, as indicated by the colored bar. The dark bar 
indicates 25 AU, for an assumed distance of 450 pc. The interferometer beam is 0.40 x 0.16 
milliarcseconds (0.18 x 0.07 AU) and the relative position uncertainties for the masers 
are smaller than the plotted symbols. We referenced the data to the strong emission at 
-3.4 kms -1 (at the origin) to stabilize the VLBA against atmospheric fluctuations along 
the line of sight. Deconvolution of the instrumental response required special attention 
because most images contained complex structure, and because of the proximity of Orion 
to the celestial equator. From our VLBA data, we estimate that the origin of the map is 
a 2000 = 5 h 35 m 14f5101 ± 0.0003, S 2000 = -05°22'30"63 ± 0.05. 

lb The distribution of H 2 masers in BN/KL within 2" of source I, observed with 
the VLA in 1983, at 0''1 resolution. The map of the SiO maser emission is shown 
within the dashed contour. The cross marks the position of source I. The registration 
uncertainty between the H 2 and SiO distributions is about 0"05, limited by uncertainty 
in the astrometric position of each with respect to the extragalactic radio reference 
frame. Arrows attached to the H 2 maser spots show the estimated relative proper 
motions over 8 years. The mean motion has been subtracted to account for motion of the 
reference emission. The double arrows indicate the direction of the proposed high-velocity 
outflow from source I. The scale bar indicates 100 AU. An arrow as long as the bar 
corresponds to a motion of about 50 kms -1 . The 1983 VLA observations covered a total 
velocity range of —12.4 to 28.1 kms -1 , except for velocities close to a 4 x 10 6 Jy H 2 
maser nare£3 with a velocity of 8 kms -1 and a position of o; 2 ooo — 5 ft- 35 m 14f 121 ± 0.003, 
^2000 = — 05 o 22'36''27 ± 0.05. The maser data were self-calibrated and relative positions 
for the maser spots measured. Absolute positions were obtained from low dynamic range 
images obtained contemporaneously by a wide-band antenna subarray that observed the 
maser flare and a phase calibrator, 0529+075. 

lc A proposed model for the region surrounding source I, showing the high- velocity and 
low-velocity outflows schematically as cones and a non-axisymmetric torus, respectively. The 
front surface of the torus is cut-away for clarity, and the SiO and H 2 maser distributions 
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are superposed. The arrows indicate schematically the high-velocity outflow; the arrows 
at position angles of —60° and 120° also correspond approximately to the direction of 
the magnetic field in OMC-l@. A toroidal geometry is suggestive but not unique. The 
H2O masers lie close to the surface of the torus, probably in outward-propagating shocks. 
Overall, the zone in which the H2O masers lie may not move significantly over time, as 
individual shocks move outward and diminish over time, while new ones are born closer to 
the protostar. Otherwise, the dynamical age of the equatorial outflow is on the order of 100 
years. 

Figure 2 A near-infrared image (3.8 /im) of the region near IRc2§, with contours 
representing 8.4 GHz emissioni and the H 2 maser distribution superposed. The bar 
indicates 1,000 AU. Labels mark the four components of IRc2 (A, B, C, and D) and other 
infrared sources. Specifically, source n is the brightest infrared object in the field. It is 
probably a young massive star with an associated HII region!, though no dust emission 
have been detected yet at 3 mm wavelength^!. The dashed box indicates the field of view 
in Figure lb. The X highlights a possible relationship bewteen the proposed outflow from 
source I and IRc2. Also shown is the center of expansion (C.O.E.) for the low-velocity 
(18 kms^ 1 ) flow in BN/KL, estimated in part from the proper motions of H 2 masers 
outside the field shownlia. 
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